Metal nanofilm in strong ultrafast optical fields 
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We predict that a metal nanofilm subjected to an ultrashort (single oscillation) optical pulse of a 
high field amplitude ~ 3 V/A at at normal incidence undergoes an ultrafast (at subcycle times < 1 fs) 
transition to a state resembling semimetal. Its reflectivity is greatly reduced, while the transmissivity 
and the optical field inside the metal are greatly increased. The deposition of the pulse energy as 
a function of field reaches its maximum and then decreases due to suppression of the active optical 
resistance. The temporal profiles of the optical fields are predicted to exhibit pronounced subcycle 
oscillations which are attributed to the Bloch oscillations and formation of the Wannier-Stark ladder 
of electronic states. These effects are promising for applications as nanoplasmonic modulators and 
field-effect transistors with petahertz bandwidth. 

PACS numbers: 71.30.+h 73.50.Fq 42. 65. Re 71.45.Gm 



Behavior of solids in strong ultrafast optical fields 
has recently attracted a great deal of attention [H-ITl). 
Such fields produce non-perturbative effects on solids, 
among which are ultrafast optical breakdown [l[ , attosec- 
ond ionization 3, 4], metallization of dielectric nanofilms 
H| > optical field-effect reversible subfemtosecond cur- 
rents in dielectrics [l(|, and electron tunneling from sur- 
faces @, 0, F° r dielectrics, when optical field is ap- 
plied with frequency fiw low enough compared with the 
band gap A vc between the valence a conduction bands, 
mostly adiabatic processes take place such as Wannicr- 
Stark (WS) localization and formation of the WS ladder 
of levels [12, HH separated by the Bloch frequency [l4| 
ujb = \e\Fa/h, where e is electron charge, F is the mag- 
nitude of the field, and a is the lattice constant. Only 
when the field F exceeds critical field F c — A.„ c /(|e|a) ~ 
2.5 V/A (for A - 10 eV and a ~ 4 A), the band gap 
is overcome by the WS splitting, and the diabatic cou- 
pling of the valence and conduction band becomes str ong , 
which can lead, in particular, to optical breakdown [lOj . 

In contrast, this Letter deals with strong optical fields 
applied to good (plasmonic) metals where there is no 
band gap at the Fermi surface and, consequently, no adi- 
abaticity for relatively low fields. In such a case, there 
are a high optical conductivity and a skin layer with a 
depth l s ~ 25 nm [l5j . Consequently for metal thickness 
h > / s , a significant or dominating fraction of the incident 
radiation energy is reflected. Interaction of the radiation 
with the metal becomes adiabatic only when the optical 
field is strong enough so > w. The plasmonic metal 
behavior seizes and WS localization [lj, [H| is estab- 
lished when, during a quarter optical period t — 7r/(2w), 
an electron acquires momentum 7r|e|-FY that exceeds the 
width 2ivh/a of the Brillouin zone. This condition is 
satisfied when the optical field F exceeds a critical field 
F c = Ahu/(\e\a) ~ 2 V/A for hu = 1.55 eV. Note that 
at this field the WS states are already strongly localized, 
Iws <^ a , where Iws = h 2 1 (ma 2 \e\F) is the WS local- 
ization radius, and m is electron mass. For F > F c , the 
strong-field regime for the metal sets on. 

As we predict in this Letter, in the strong-field regime 



the optical properties of the metal differ dramatically 
from those at low to moderate fields, becoming remi- 
niscent of scmimetals: plasmonic properties and strong 
reflection associated with the skin effect are suppressed 
during subcycle time intervals driven by the instanta- 
neous optical field. Light transmission through the metal 
is increased but the optical absorption in the metal is re- 
duced at very high fields. Ultrafast behavior of the metal 
is radically changed: both the reflection and transmis- 
sion exhibit subcycle Bloch-type oscillations with period 
t f ~ 2-k/cjb; e.g., t f ~ 0.5 fs for F = 2.5 V/A. Thus, 
the strong optical field is predicted to cause ultrafast re- 
sponse of the metal irrespectively of the frequency and 
duration of the excitation pulse. Under these conditions, 
the Bloch oscillations are predicted to manifest them- 
selves in natural metals, while earlier such oscillations 
were observed only in artificial semiconductor superlat- 
tices - see, e.g., Refs.[l6l-[T8i These predicted strong- field 
effects open up routes toward using metals as active ele- 
ments for deep ultrafast modulation of optical fields - cf. 
the perturbative modulation in active plasmonics [19(. 

Consider an ultrashort optical pulse incident normally 
on a metal nanofilm. Propagation of such a pulse is de- 
scribed by the Maxwell equations where dielectric polar- 
ization is determined by quantum dynamics of the elec- 
trons in the metal. This dynamics is described by the 
Schrodinger equation in the presence of the electric field 
inside the metal. We neglect the Coulomb scattering of 
electrons because the characteristic time r s of such scat- 
tering in metals significantly exceeds the l eng th of our 
optical pulse (e.g., t s = 20 — 50 fs in silver [la]). 

We solved numerically the coupled system of the 
Maxwell and Schrodinger equations usin g th e finite dif- 
ference time domain (FDTD) method [20l l2l| for a finite- 
size system with the absorbing boundary conditions for 
Maxwell equations. The tight-binding model was used in 
the solution of the Scrodinger equation. The size of the 
computational space in the direction of propagation of 
the pulse [z direction) was 6000 nm. The metal film was 
placed at the midplane of the system, i.e., at z = 0. In 
numerical solution of the Maxwell equations, we assumed 
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that the spatial step was 1 nm and the time step was 0.7 
attoseconds (1 as=10~ 18 s). The optical pulse was gen- 
erated at the left boundary and propagated along the 
positive direction of the z axis with the polarization of 
the electric field along the x axis. 

We assume that a single-oscillation pulse form, 



F x {t) = F e- u2 (l-2u 2 



(1) 



where Fq is the amplitude, which is related to power 
V = cFq/Att, where c is speed of light, u = t/r, and 
r is the pulse length, which is set r = 1 fs in our calcula- 
tions. Similarly short pulses have been recently used [Toj . 
The metal can be described by one-particle Schrodinger 
equation with the Hamiltonian 
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V(r) + eF x (z,t)x, 



(2) 



where V(r) is the periodic crystal potential, and F x (z,t) 
is the electric field inside the metal, which is found from 
the solution of the Maxwell equations. Without the elec- 
tric field, the electron system has standard band struc- 
ture. Below we consider one conduction band (CB, or 
sp-band in silver) and one valence band (VB, or d-band). 

We assume that the periodic potential V(r) is sepa- 
rable in all three directions, x, y, and z, with period a. 
Then the electron dynamics along the direction x of elec- 
tric field separates. For each ba nd, the energy dispersion 
law has the tight-binding form (22|, HH 



E a {k) — e a H — ^ cos(fca), 



(3) 



where a — c or v for CB and VB, respectively, A a is the 
width of band a, and e a is the band offset. In the absence 
of the optical field, the wave functions satisfy the Bloch 
theorem. 



1 



ipak{x) = —e l x u ak (x) , 

Z7T 



(4) 



where u ak {x + a) — u ak {z) are periodic Bloch unit-cell 
functions, and k is the (pseudo) wave vector. 

In the presence of the optical field, F x (z, t), we express 
the general solution of the time-dependent Schrodinger 
equation in the basis of the Bloch functions as 



9(x,z,t) 



£ 



7r /a 



dk <f>a{k,Z,t)lp a k(x) , (5) 



— 7r/ a 



Here the time-dependent wave vector is defined as 



k T (q,t) 



F x (z,ti)dti 



(8) 



and S(k) — J2 n ^(^ + 2im/a), where n — 0, ±1, . . . , and 
5(k) is the Dirac delta- function. The Houston functions 
are exact solutions of the time-dependent Schrodinger 
equation for a single band with the Bloch function ip aq (x) 
((D) as the initial condition at t — — oo. 

Expansion coefficients /3 a (q,z,t) satisfy the equations 



d/3 a (q,z,t) 
dt 



F x (z,t) 



^2 Q aa >(q,z,t)(3 a ,(q,z,t) 



(9) 



where we denoted 

Qaa'(q,Z,t) = Z ao - exp <j / 

A a ~A a : r ' 
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dz u ak (z)*i—u a , k (z) 



F x (z,t 2 )dt 2 



(10) 



Here dipole matrix elements Z aa i describe diabatic cou- 
pling of VB and CB in optical field. 

The electric current generated by electron dynamics 
in the optical electric field pulse has two contributions: 
interband and intraband, is J x = J x ntcl + J> ntra . The 
interband polarization has the following form 



P™ ter (z,t) 
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2ira 3 



dq E fM) \B M Hq,z,t)Q( q ,z,t)B^(q,z,t) 



fl—V,C 



(11) 



where Q is matrix with elements (JTUJ) Q aa i , and = 
($ v ,$ c ) is a two-component vector, which is determined 
by the solution of Eq. © with the following initial con- 
ditions: = (1,0) and B^ = (0,1). Here f^q) is 
the Fermi factor, which is 1 for initially occupied states, 
i-e. fn=v{q) = 1 and f^ =c {\q\ < k F ) = 1, where k F 
is the Fermi wave vector, and it is zero otherwise. The 
interband current is 4 ntor (z,£) = dP x Dtor (z,t)/dt. The 
intraband current is due to shifting of electrons and is 



where <fi a (k, z, t) can be expressed in term of the Houston 
functions HU ^Jp(k,z,i), 

Mk,z,t) = ^i3 a (.q,z,t)$Vp(k,z,t) , (6) 
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$i H q Xk,z,t) = ~5(k-k T (q,t)) x 

cxp j-i (tjr- + ~^f-J dh cos[kT(q,ti)a]\^ . (7) 



Jt^{z,t) 



£ U(q) 



2na 3 



dq sin [fcr(<7, t)a\ x 



2h 



(12) 



Below we use model parameters corresponding to the 
band structure of silver: e v = —4.1 eV, e c = eV, = 



3 
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FIG. 1: Reflected and transmitted pulses, (a) Spatial distri- 
butions of the electric field as functions of the propagation co- 
ordinate z shown for different values of Fq. The metal film of 
thickness 25 nm is placed at the center [z = 0) and depicted 
as the red stripe. The distribution of electric field consists 
of the reflected (to the left) and transmitted (to the right) 
pulses propagating in the opposite directions. The size of the 
computational field in the z direction is 6000 nm. (b) Spec- 
tral intensities of reflected, transmitted, and sum pulses for 
Fo = 4.5 V/A as functions of optical frequency / = uj/(2ty). 
Fragments of the curves at / > 2 PHz are also shown with 
the x 10 magnification as indicated on the graph. 



0.82 eV, and A c = -9.1 cV. We choose Z vc = 0.7 eA 
[2o| . The thickness of the film is set ft. = 25 nm. 

In Fig. Q] (a) , the spatial distribution of the pulse elec- 
tric field is shown for the reflected (left) and transmitted 
(right) pulses for amplitudes Fq. For a relatively small 
field, Fo = 1 V/A < F c , the nanofilm behaves as a reg- 
ular metal with a pronounced skin effect and strong re- 
flection of the incident pulse. With increasing the field, 
F = 2.5 V/A > F c , the response of the electron system 
is highly nonlinear, and the metal film becomes relatively 
transparent. Both the reflected and transmitted pulses 
are strongly reshaped compared to the incident pulse. As 
the pulse peak field further increased to Fq = 3.5 V/A 
and Fo = 4.5 V/A, the film transparency is further in- 
creased, in a sharp contrast to metallic behavior. Also, 
there are pronounced sub-wavelength oscillations in the 
pulse shape for both the reflected and transmitted fields, 
which we interpret as an effect of the Bloch oscillations. 

An important peculiarity of Fig. [T] (a) is a nonzero 
area of both the transmitted and reflected pulses, 9 f ' r — 
J^ oo F(z t ' r ,t)dt ^ 0, where z* > and z r < 0, and 
t and r denote transmission and reflection, correspond- 
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FIG. 2: The electric field of the incident pulse (black lines) 
and the electric field at the midpoint of the metal film (red 
line) are shown for different values of the peak electric field 
Fo of the incident pulse: (a) F = 2.0 V/A, (b) F = 2.5 V/A, 
and (c) F = 3.5 V/A. 



ingly. This imbalance of the field integrals is due to 
high nonlinearity of the field interaction with the metal. 
When the absorption in the matter is small, which is the 
case presently, then the total area still is "near zero", 
|0 r + 0*| < \6 r \ + Note that pulses with 6 ^ 

were considered previously in nonlinear optics [26[. For 
ultrashort pulses, the magnitude and sign of Q t,r is deter- 
mined by the carrier-envelope phase as characteristic for 
any nonlinear effects induced by few-oscillation pulses, 
cf. Refs. llCj, |l_l|. (For the pulses in the present Letter, 
the carrier-envelope phase is zero.) 

Other important features of the transmitted and re- 
flected pulses are illustrated in Fig. [1] (b) displaying 

2 

their spectral intensities /y r = F^ r , where index 

/ = uj / (2tt) indicates Fourier transform in terms of lin- 
ear frequency. Besides peaks at the carrier frequency, 
/ w 0.25 PHz, there are also pronounced peaks at the 
"third harmonic" frequency, / « 0.75 PHz, that are due 
to the strong nonlinearity of the optical response at such 
fields. There are also appreciable peaks at the Bloch fre- 
quency, / ~ 3.5 PHz, which are shown magnified by a 
factor of xlO. Observation of these peaks would be the 
first evidence of the Bloch oscillations in real crystals. 

This enhanced transmission of the ultrastrong optical 
pulse is accompanied by an increase of the electric field 
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Peak electric field, F (V/A) 



FIG. 3: (a) The reflectance of optical pulse (black line) and 
the maximum electric field at the midpoint of the nanofilm 
(red line) are shown as functions of the peak electric field Fo 
of the incident pulse. The maximum electric field in the metal 
film is shown in units of the peak electric field Fo. (b) The 
absorbance of the optical pulse is shown as a function of the 
peak electric field Fq. 
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FIG. 4: Bloch oscillations in the transmitted field and elec- 
tron momentum, (a) Electric field distribution in space of 
the transmitted optical pulse is shown for different values 
of Fo. (b) The corresponding dimensionless time-dependent 
wave vectors in the first Brillouin zone akr{q = 0, t) as func- 
tions of time t. The origin of time is chosen arbitrary, and 
the graphs are offset vertically for clarity. 



F metal inside the metal nanofilm. In Fig. [2 we show 
the time evolution of this field at the midplane of the 
nanofilm in comparison to the incident pulse. At a near- 
critical pulse field F = 2 V/A ~ F c , see Fig. H (a), the 
electric field F mGtal inside the metal is much weaker than 
that of the incident pulse. In contrast, for larger pulse 
amplitudes (F = 2.5,3.5 V/A), see Fig. M 0, c), the 
electric field F mctal becomes comparable to the incident- 
pulse electric field. The sharp peaks and high-frequency 
oscillations of electric field F metal are due to the Bloch 
oscillations in electron dynamics, as we discuss below. 

The reflectance of the optical pulse (a fraction of the 
reflected pulse energy) is shown in Fig. [3]Ja) as a func- 
tion of Fo. Strong suppression of the pulse reflectance 
and correspondingly increase of the pulse transmission 
for F > F c ~ 2 V/A is clearly visible. The suppression 
of the reflectance is correlated with increase of the electric 
field inside nanofilm shown by the red line in Fig. [3^a) as 
F™^ al /F . The strong chan ges in both reflectance and 
internal electric field occur at Fo ~ F c ~ 2 V/A. 

The absorbance of the metal nanofilm, calculated 
as a fraction of the pulse energy dissipated inside the 
nanofilm, is illustrated in Fig. [3] (b) as a function of 
amplitude Fo. This predicted behavior is very unusual. 
At a low pulse amplitude, the absorbance is under- 
standably low due to the skin effect since most of the 
pulse is reflected back. Then as Fo increases, the ab- 
sorbance increases dramatically reaching w 4% maximum 
for F « 3 V/A > F c , which is attributed to the Bloch 
oscillations and WS localization leading to the collapse 
of the skin effect - cf. the red line in panel (a). With 
the further increase of Fo, the absorbance decreases de- 
spite the field in the metal being almost the same. This 
can be understood as the nanofilm acquiring properties 
of semimetal with a low active conductivity and, conse- 
quently, low optical loss characteristic of semimetals in 
contrast to good (plasmonic) metals. 

The power density dissipated from a single pulse in the 
metal reaches its maximum also at Fo « 3 V/A. After its 
dissipation and thermal equilibration, this causes an esti- 
mated increase of the nanofilm temperature by ~ 500 K. 
Thus the metal may survive such a high-intensity pulse 
without a damage. This is attributed to loss of the metal- 
lic and onset of semimetallic properties by the nanofilm. 
Our theory explains the absence of optical damage of 
metal electrodes subjected to comparable pulses in ex- 
periments [1(| . A resistance of dielectrics to high power 
intensities has also been proven experimentally [l(| ■ 

The origin of this highly nonlinear behavior of a metal 
film in a strong optical field can be understood from elec- 
tron dynamics within a single conduction band. In the 
optical field, an electron with initial wave vector q is mov- 
ing in the reciprocal space following the time-dependent 
wave vector kT(q,t), see Eq. ©. Therefore, all electrons 
are shifted in the reciprocal space by the same wave vec- 
tor Aq(t) = | J* F x [z,t\)dt\, and the net current is gen- 
erated. For a strong pulse, the shift Aq is large and, 
for Fq > F c , can become greater than the Brillouin zone 
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extension k = 2ir/a, causing the Bragg reflection of the 
electrons, which results in the Bloch oscillations and in 
the onset of the WS localization states. 

^ Fq ^ Fci the electron current acquires oscillations 
at Bloch frequency u>b > 4w, which suppresses the sus- 
ceptibility at the optical frequency cj. This results in the 
loss of the metallic optical properties. Since the electron 
spectrum is discreet (the WS ladder) and, consequently, 
the density of states at the Fermi level is zero, the metal 
in the strong optical fields behaves as a semimetal with 
a relatively high transparency and low reflection. 

In Fig. |U[a) , the spatial distribution of electric field in 
the transmitted optical pulse is shown for different val- 
ues amplitude Fq. With increasing Fq above the thresh- 
old F c w 2 V/A, well pronounced Bloch oscillations de- 
velop in the field distribution. Their total number is pro- 
portional to field amplitude, n w \e\aFo / (2huj) . This is 
the number of times that an accelerated electron crosses 
the Brillouin zone boundary, as can be illustrated by the 
comparison to the temporal dependence of the electron 
quasi- momentum displayed in Fig. SJb). These Bloch 
oscillations are also visible in the temporal evolution of 
electric field inside the metal film - Fig. [5J 

To briefly conclude, we have predicted a highly un- 
usual and interesting behavior of metal nanofilms sub- 



jected to strong ultrashort (single-oscillation) optical 
pulses with the field amplitude ~ 3 V/A (intensity 
~ 2.4 x I0 14 W/cm 2 ). This includes such effects as 
reduction of the metallic high reflection (suppression of 
the skin effect) and greatly increased transmission of the 
pulse energy through the nanofilm, while both the ab- 
sorbance and the total energy deposition dramatically de- 
crease at the high pulse intensity. This indicates that the 
optical field induced a transition to a semimetallic state. 
These phenomena develop at subcycle times < 1 fs and 
are driven by the pulse instantaneous amplitude. The 
pulse shapes of the reflected, transmitted, and internal 
fields hint that the metal is likely to return to its original 
state by the end of the pulse. Thus this transformation of 
the metal to a semimetal is optical-field-induced and re- 
versible. The modulation depth for the transmitted pulse 
is extremely high, which shows prospects of using metals 
as active elements in ultrafast modulators and field-effect 
transistors with petahertz bandwidth. 
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